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ABSTRACT 

I present a predictive analysis for the behavior of the far-infrared (FIR)-radio correlation as a 
function of redshift in light of the deep radio continuum surveys which may become possible using 
the square kilometer array (SKA). To keep a fixed ratio between the FIR and predominantly non- 
thermal radio continuum emission of a normal star-forming galaxy, whose cosmic-ray (CR) electrons 
typically lose most of their energy to synchrotron radiation and Inverse Compton (IC) scattering, 
requires a nearly constant ratio between galaxy magnetic field and radiation field energy densities. 
While the additional term of IC losses off of the cosmic microwave background (CMB) is negligible 
in the local Universe, the rapid increase in the strength of the CMB energy density (i.e. ~ (1 + z) A ) 
suggests that evolution in the FIR-radio correlation should occur with infrared (IR; 8— 1000 (im) /radio 
ratios increasing with redshift. This signature should be especially apparent once beyond z ~ 3 
where the magnetic field of a normal star-forming galaxy must be ~50 fiG to save the FIR-radio 
correlation. At present, observations do not show such a trend with redshift; z ~ 6 radio-quiet quasars 
appear to lie on the local FIR-radio correlation while a sample of z <~ 4.4 and z ~ 2.2 submillimctcr 
galaxies (SMGs) exhibit ratios that are a factor of ~2.5 below the canonical value. I also derive a 
5cr point-source sensitivity goal of ~20 nJy (i.e. ctrms ~ 4 nJy) requiring that the SKA specified be 
A e s/T sys s» 15000 m 2 K _1 ; achieving this sensitivity should enable the detection of galaxies forming 
stars at a rate of > 25 M Q yr _1 , such as typical luminous infrared galaxies (i.e. Lir > 10 11 L Q ), at 
all redshifts if present. By taking advantage of the fact that the non-thermal component of a galaxy's 
radio continuum emission will be quickly suppressed by IC losses off of the CMB, leaving only the 
thermal (free-free) component, I argue that deep radio continuum surveys at frequencies >10 GHz 
may prove to be the best probe for characterizing the high-z star formation history of the Universe 
unbiased by dust. 

Subject headings: galaxies: evolution - radio continuum: galaxies - magnetic fields - infrared: galaxies 



1. INTRODUCTION 

Radio continuum emission from galaxies arises due to 
a combination of thermal and non-thermal processes pri- 
marily associated with the birth and death of young mas- 
sive stars, respectively. The thermal (free-free) radiation 
of a star-forming galaxy is emitted from Hn regions and is 
directly proportional to the photoionization rate of young 
massive stars. Since emission at GHz frequencies is op- 
tically thin, the thermal radio continuum emission from 
galaxies is a very good diagnostic of a galaxy's massive 
star formation rate. 

Massive (> 8 Mq) stars which dominate the Lyman 
continuum luminosity also end their lives as supernovae 
(SNe) whose remnants (SNRs) are responsible for the ac- 
celeration of cosmic-ray (CR) electrons into a galaxy's 
general magnetic field resulting in diffuse synchrotron 
emission. Thus, in a more complicated manner, the non- 
thermal radio continuum emission also traces the most 
recent star formation activity in a galaxy. Ideally, one 
would like to isolate the thermal component as it is a 
more direct measure of the most recent massive star for- 
mation activity. However, at GHz frequencies, the non- 
thermal fraction typically dominates the total radio con- 
tinuu m emission (i.e. ~10:1 at ~1 GHz: [Condon fc Yin I 
1990) making the isolation of the thermal fraction diffi- 
cult. 

These same massive stars are often the primary 



sources of dust heating in the interstellar medium 
(ISM) as their starlight is absorbed and reradiated 
at far-infrared (FIR) wavelengths by interstellar 
grains. This common origin between the FIR dust 
emission and thermal + non-thermal radio con- 
tinuum emission from galaxies is thought to be 
the dominant physical processe s driving the FIR- 



radio correlation on global (e.g. Ide Jong et al. 1 1 1985 
Helou. Soifer. fc Rowan-Robinson 1985: iNiklas I 11997 : 
Niklas fc Bec k I 1 19971: lYun. Reddv. fc Condon I 12001 ) 
and local ( e.g. Beck fc Golla I Il988t IXu et al. I I199S : 
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2008; iHughes et aTT feooe) scales. While massive star 
formation provides a shared origin between emission at 
these two wavelengths, a number of physical processes 
must conspire to yield such a remarkably constant ratio 
among galaxies s panning nearly 5 orders of mag nitude in 
luminosity (e.g. lYun. Reddv. fc Condonl 120011 ). among 
other properties (e.g. Hubble type, FIR color, and 
FIR/optical ratio). For example, since the non-thermal 
synchrotron emissivity of a single electron is roughly 
proportional to the electron density times the square 
of the magnetic field strength, it would seem that the 
FIR-radio correlation should be highly dependent on the 
propagation of CR electrons and a galaxy's magnetic field 
distribution and strength. If magnetic fields in galaxies 
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are built-up over time, and therefore were weaker in 
the past, the correlation should be different at higher 
redshifts. However, all indications suggest that the FIR- 
radio correlation holds out to moderate redshifts (e.g. 



Garrett I I2002t iGruppioni et al. I [20031: lAppleton et al. 



20041: iFraver ■ et al.T 120061: iMurphv et al 



2009a: 



M.T Sargent et al. I 12009. in preparation! ) suggesting 



that the magnetic field strength and structure in these 
galaxies is similar to what is observed for galaxies 
in the local Universe. Consequently, comparing the 
FIR and radio emission characteristics of galaxies can 
illuminate properties that are typically inaccessible by 
other observational means. 

At present, detailed studies of the FIR-radio corre- 
lation among star-forming galaxies has been limited to 
only moderate redshifts at z < 1. With the sensitivity 
of existing FIR and radio capabilities, studies at higher 
redshifts only probe the most extreme objects. For ex- 
ample, deep (o~ rms ~ 0.55 mJy) 70/im observations 
(|Fraver et al. I 2006) from the the Far-Infrared Deep 
Extragalactic Legacy (FIDEL; PI. M. Dickinson) sur- 
vey is able to probe luminous infrared galaxies (LIRGs; 
10 11 < Lir < 10 12 Lq) and ultraluminous infrared 
galaxies (ULIRGs; > 10 12 L Q ) out to redshifts of z ~ 1 
and z ~ 2, respectively. With Herschel, surveys such as 
GOODS- Her schel should improve this situation by mea- 
suring the FIR properties of normal star-forming galax- 
ies out to z ~ 1, and that of LIRGs and ULIRGs out to 
redshifts of z ~ 2 and z ~ 4, respectively. ALMA will 
also play a significant role by measuring the peak of the 
FIR spectral energy distribution (SED) for all LIRGs at 
z > 5. 

These populations of dusty star-forming galaxies ap- 
pear to dominate the stellar mass assembly at in- 
creasing redshifts; the star formation rate density in- 
creases by a factor of ~ 5—10 between z ~ 
and z ~ 1, becoming increasingly obscured (i.e. > 
60%) by dust (lElbaz et al. H2002t ICharv fc Elbaz l l200lt 
iLe Floc'h et al. 1 120051 : iMagnelh et. al. I I2009D . Thus. 
LIRGs and ULIRGs appear to dominate the luminos- 
ity density at increasing redshifts, and optically thin 
measures of star formation and active galactic nuclei 
(AGN) activity are critical for proper quantification of 
stellar mass build-up over cosmic time. At radio wave- 
lengths, however, detecting such high redshift galax- 
ies remains extremely difficult. Even with a fully op- 
erational EVLA, IR-bright star-forming galaxies (e.g. 
M 82; Lib, m 4 x 10 10 L ) and moderate LIRGs (i.e. 
L IR w 3 x 10 11 Lq) will not be detectable beyond red- 
shifts of z ~ 1 and z ~ 2, respectively. 

A next-generation radio facility such as the Square 
Kilometer Array (SKA) should easily remedy this dis- 
parity between the depth of FIR and radio continuum 
surveys. While the SKA was initially proposed solely 
on the basis of Hi science, deep continuum imaging is 
critical for the realization of nearly all of the five estab- 
lished Key Science Projects (KSPs), especially studies of 
"The Origin and Evolution of Cosmic Magnetism" and 
"Galaxy Evolution and Cosmology." Essential to these 
two science goals is the proper measurement of the star 
formation and AGN history over cosmic time for which 
deep radio continuum studies may provide an excellent 
advantage over other wavelengths. 



In this paper I present physically motivated expecta- 
tions for the behavior of the FIR-radio correlation at in- 
creasing redshift and discuss their implications. I also 
discuss how the SKA, when combined with future FIR 
observations to be obtained with next generation facil- 
ities both prior to, and commensurate with, final SKA 
science operations, will be able to tackle interesting prob- 
lems associated with these KSPs. The paper is orga- 
nized as follows: In §2 I discuss the FIR-radio correla- 
tion and introduce the physical processes for which it 
depends on. Then, in section §3, I present the results 
for the expected evolution in the FIR-radio correlation 
at increasingly high redshifts and report on the technical 
specifications for detecting high-z galaxies in deep radio 
continuum surveys. In Section §4 I discuss the physical 
implications of the results for characterizing the proper- 
ties of galaxies at high z and compare these theoretical 
expectations with existing observations. Finally, in §5, I 
summarize my conclusions. 

2. THE FIR-RADIO CORRELATION: PHYSICAL 
CONSIDERATIONS 

A m ajor result of the Infrare d Astronomical Satellite 
(IRAS: iNeugebauer. et aTlll984h all-sky survey was the 
discovery of a correlation between the globally measured 
far infrared (FIR; 42 — 122 fim) dust emission and the 
optically thin radio continuum emission of normal late- 
type s tar-forming galaxies without AGN (de Jong et al. I 



[19851 : IHelou. Softer, fc Rowan-Robinson I I1985T T The 
most remarkable feature of this correlation is that it 
displays such little scatter (i.e. «0.26 dex) among 
galaxies spanning 5 orders of magnitude in luminosity 
(lYun. Reddv. & Condon I l200ll) . While the FIR emis- 
sion is due to the thermal re-radiation of interstellar 
starlight by dust grains, the radio emission is primar- 
ily non-thermal synchrotron emission from CR electrons 
that propagate in a galaxy's magnetic field after initially 
being accelerated by SN shocks or other processes. The 
physics that maintains a strong correlation between these 
two quantities over such a wide range of galaxies remains 
unclear, however departures from the nominal ratio can 
shed insight on a number of typically inaccessible galaxy 
properties. 

While the FIR emission is dependent upon the chem- 
ical make-up and size distribution of dust grains, which 
is hard to decipher, a galaxy's non-thermal radio emis- 
sion is relatively simple to interpret. What is not sim- 
ple, however, is determining which energy-loss processes 
may dominate a galaxy's population of CR electrons, as 
well as the role of escape. CR electrons will cool via 
a number of mechanisms as they propagate through the 
ISM of galaxies. These energy-loss processes include syn- 
chrotron radiation and Inverse Compton (IC) scattering, 
as well as ionization, bremsstrahlung, and adiabatic ex- 
pansion losses. Below I briefly describe the role that each 
of these energy-loss processes may play for the cases of 
norma l star-forming ga laxies and compact starbursts fol- 
lowing [Longair] (|1994h . 

2.1. Essential Physics: Normal Star-Forming Galaxies 

For a typical star- forming galaxy, CR electrons primar- 
ily lo se their energy due to synchrotron and IC processes 
(e.g. lCondon~|[l99a ). although escape also plays a role 
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Fig. 1. — Left: The cooling timcscales of CR electrons as a function of electron energy due to a number of physical processes: synchrotron, 
Inverse Compton (IC), ionization, and bremsstrahlung. For these calculations typical parameters for normal star- forming galaxies are 
assumed: B = 10 fiQ, f7 ra d = 10 -12 erg cm" 3 , and nisM = 1 cm" 3 . A timescale due to the escape of CR electrons from a galaxy, assuming 
random walk diffusion and an escape scale-length of l csc = 3 kpc, is also included. The top axis indicates the synchrotron emitting frequency 
of the CR electrons. The shaded region indicates the electron energies emitting predominantly in the 0.5 — 2 GHz passband. Right: Same 
as what is plotted in the left panel, but for values that may be more typical of galaxies hosting compact starbursts, namely: B = 1 mG 
and njgM = 10 4 cm -3 , which assumes the usual B oc ^/niSM scaling. The additional term from adiabatic expansion losses due to a 
starburst-driven galactic wind, having a velocity of v w = 300 km s" 1 and assuming a disk scale-height h = 1 kpc, are also shown. L/ ra d was 
also scaled up to 10 — 8 erg cm -3 . For normal star-forming galaxies, synchrotron and IC losses dominate CR electron energy losses, along 
with escape, for radiation observed at >GHz frequencies. However, in the case where a mG magnetic field is considered, synchrotron and 

IC losses no longer are the dominant cooling processes for CR electrons emitting at GHz frequencies. 
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the magnetic field energy density of the galaxy. Using 
Equation [TJ the synchrotron cooling timescale, 
E /\dE /dt\ syn , for CR electrons can be expressed as 

1/5! / t-? \ ^ / 



' syn 

yr 



-1/2 
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syn 



Up 
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12 erg cm 3 



(2) 

Naturally, the synchrotron cooling timescale for CR elec- 
trons can be rewritten as 
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Fig. 2. — A model radio-FIR spectrum for a galaxy having a FIR 
flux of FpiR = 10" 13 W m~ 2 . The FIR portion of the spectrum is 
given using the IR dust emission component of the I Dale & Helou I 
(2002) SED libraries. The thermal radio spectrum, which goes 
as Sj oc v~ 01 , is normalized assuming that the ionizing photon 
rate estimated by the thermal radio continuum and IR luminosity 
are the same, and that the thermal fraction at 1.4 GHz is 10%. 
Under these assumptions the nominal FIR-radio correlation (i.e. 
gjR ^ 2.64) is achieved. The non-thermal spectrum includes energy 
losses from all of those mechanisms considered in deriving the CR 
cooling timescales shown Figurefl] (See §2.1 for details). 



Similarly, the energy loss of CR electrons due to IC scat- 
tering goes as dE/dt oc U ra aE 2 where U Ia ^ is the radia- 
tion field energy density of the galaxy. Equation [1] can 
therefore be used again to write the IC cooling timescale, 
tig = E/\dE/dt\ IC , as 
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(e.g. lHelou fc Bicav lfl993T ). Let us assume that CR elec- 
trons propagating with a pitch angle <f> in a magnetic 
field of strength B have isotropically distributed veloci- 
ties such that < sin 2 4> >= § leading to B± w 0.82 B. A 
CR electron having energy E will emit most of its energy 
at a critical frequency v c where 
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For GeV electrons considered here, the bulk of losses arise 
from interactions with IR/optical photons, which domi- 
nate [Trad- 

The effective cooling timescale for CR electrons due to 
synchrotron and IC losses is 

T cool = T syn + T IO ' (^) 

which, by combining Equations [5] and 01 we can express 
as 

7 / V, \ -1/2 / R\'/ 2 

; 5.7xl0 7 ( 
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The energy loss of CR electrons by synchrotron radi- 
ation goes as dE/dt oc UbE 2 , where Ub — B 2 /(8ir) is 



(6) 

If CR electrons do not lose all of their energy as they 
propagate through the ISM of a galaxy, they may even- 



4 



MURPHY 








6 



8 



10 



12 



Fig. 3. — The expected observed-frame 1.4 GHz flux density for galaxies of varying IR luminosities assuming the FIR-radio correlation 
(i.e. qiR, = 2.64) as a function of redshift for nominal galaxy parameters (Table [TJ. Additional energy losses to CR electrons arising from 
IC scattering off of the CMB, whose radiation energy density scales as ~ (1 + z) 4 , are included. Flux densities, corrected for bandwidth 
compression, are plotted assuming internal magnetic field strengths of 10 fiG (solid-line) and 100 fiG (dotted-line); at high z the non-thermal 
emission from galaxies is highly suppressed even if B ~ 100 /iG, leaving only the thermal component detectable. The case of a moderate 
LIRG (iiR ~ 3 X 10 11 Lq, implying a star formation rate of ~ 50 Mq yr _1 ) is highlighted for which a 50 fiG magnetic field strength 
(dashed-line) is also shown. The sensitivities of existing deep radio continuum surveys using the VLA, the expected depth of the EVLA, 
and suggested depth for the SKA (Table |2}, which would detect typical LIRGs (i.e. galaxies forming stars at rates > 25 Mq yr" 1 ) at all 
z if present, are shown. At this depth the SKA will be sensitive to galaxies lying in the shaded region which indicates the e xpected flux 
density range for galaxies included in the z ~ 4, 5, 6, 7, and 9 UV luminosity function studies of Bouwcns ct al. (2007, 2008). 



tually escape the system. In simple diffusion models, 
the propagation of CR electrons is usually character- 
ized by an empi rical, energy-dependen t diffusion coef- 
ficient, De (e.g. iGinzburg et al. Ill980f ). Values of De 
has been found to be around 4 — 6 x 10 28 cm 2 s _1 for 
<GeV CRs by fitting diffusion models with direct mea- 
surements of CR nuclei (i.e. secondary-to-primary ra- 
tios like B oron-to-Carbon) w i thin the Solar Neighbor- 
hood (e.g. Pones et al. 1120011: iMoskalenko et al. 1 120021 : 



iMaurin. Taillet fc Donato I I2002P )" While this empiri- 
cally measured value is for CR nuclei within the Milky 
Way, it has been found to be consistent with inferred dif- 
fusion coefficients for CR electrons both radially along 
field lines (~ 10 29 cm 2 s _1 ) and vertically across field 
lines (~ 10 28 cm 2 s _1 ) in t he thin disks of galaxies 



NGC 89 1 and NGC 4631: [Dahlem. Lisenfeld. fc Gollal 
119951) . iHeesen et al~l ([20090 also report a diffusion 
coefficient of ~ 2 x 10 29 cm 2 s _1 in the halo of 
NGC 253. Similar values near ~ 10 29 cm 2 s _1 are 
also found in hydrodynamic simulations of bubble/ super- 
bubb l e induced galaxy ou tflows (e.g. iRebusco et al. I 
120051: iRoediger et al. I 12007m . Using the average value 
of the diffusion coefficient obtained through direct mea- 



surements of CR nuclei, which is quite similar to that 
inferred for radial and vertical diffusion coefficients for 
CR electrons, the diffusion coefficient is expressed as 
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5 x 10 28 , 
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E < 1 GeV 
i 1 / 2 , E > 1 GeV. 



(7) 



Then, for for the case of random walk diffusion, such that 
Tcsc = le SC /DE, where l esc is the escape scale-length, CR 
electrons will escape a galaxy in 
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, E > 1 GeV, 



where Equation [T] has been used to express the electron 
energy dependence in terms of its synchrotron emitting 
frequency for a given magnetic field strength. 

Most edge-on galaxies possess thick synchrotron 
disks with scale- heights ranging between 1 — 3 kpc 
(|Lisenfeld. Dahlem. fc Ehlel I2004D . In the case of 
equipartition between the energy density of cosmic rays 
and the magnetic field, the CR scale-height will be 
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(3 + a)/ 2 w 2 times that of the synchrotron scale- 
height for a synchrotron spectral in dex of a w 0.8 (e.g. 
INiklas. Klein, fc WielebinskTl [l997h . Then, assuming a 
magnetic field strength of 10 fiG and a CR scale-height 
of ~3 kpc (synchrotron scale-height of 1.5 kpc), the typ- 
ical escape time of a CR electron emitting at 1.4 GHz is 
~3 x 10 7 yr. This is comparable to the expected cooling 
time of such electrons to synchrotron and IC processes, 
assuming that U la d — 10 -12 erg cm -3 , indicating that 
normal galaxies do not behave like calorimeters. 

Vertical diffusion of CRs, eventually leading to their 
escape into intergalactic space due t o either open mag- 
netic field lines or Par ker instabilities (lParke"r"lll966l ) , was 
looked at in detail bv lHelou fc Bicavl ( 19931 ). Their pro- 
posed physical model for the FIR-radio correlation in- 
cludes a leaky-box confinement scheme in which the es- 
cape scale-length of CR electrons is nearly independent 
of ISM density and magnetic field strength, but scales 
with disk scale-height to keep a constant IR/radio ratio. 
While escape may play a role, injecting scatter into the 
observed IR/radio ratios in larger spirals, it is notewor- 
thy that for the case of irregular galaxies, which generally 
lack dense ISM and magnetic field to keep CR electrons 
trapped once leaving their initial clouds around SNRs, 
escape appears to play a much more criti cal role pushing 
such galax i es off of the correlat ion (e.g. iCannon et al. I 
[20051 120061 : fMurphv et al. Il200l . 

2.2. Essential Physics: Starbursting Galaxies 

Additional energy-loss terms that may become increas- 
ingly important in the case of galaxies hosting strong 
starbursts are now considered. In such systems, whose 
energetics and ISM may be vastly different, ionization, 
bremsstrahlung, and adiabatic cooling through advection 
out of a galaxy by galactic-scale winds, can all play a 
significant role in cooling CR electrons. Ionization losses 
begin to dominate over synchrotron and IC losses for CR 
electron energies below about ~1.3 GeV such that, 
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Combining this with Equation [T] yields, 
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Thus, if magnetic fields are sufficiently large in 
sta rbursting systems (i.e. >1 mG), as suggested 
by iThompson et al.~l (|2006f ). electrons radiating GHz 
synchrotron continuum emission will have energies 
<0.3 GeV and ionization losses should dominate over 
synchrotron cooling likely leading to depressed syn- 
chrotron emission (see Figure 1). 

Furthermore, bremsstrahlung losses are probably negli- 
gible in normal galaxies given that the CR elec tron mean 
free path is ~5 g cm -2 ([Garcia-Munoz |[l977t) . an order 
of magnitude smaller than the ~60 g cm -2 for the radia- 
tion length of the ISM. However, for a sufficiently dense 
ISM, bremsstrahlung losses may exceed ionization losses 
for electrons at energies above ~1.1 GeV such that, for 



(li^ ^ 8.6 x 10 7 P^f 1 . (11) 
\ yr J Vcm-- 5 / 

Additionally, if these starbursts are compact, the adia- 
batic expansion losses to CR electrons may become more 
rapid due to the advection out of each system with a 
galactic wind, having speed v Wl where the adiabatic life- 
time goes as, 



Tad 



1.0 x 10 a 



h 
kpc 



Vkm s- 1 / 



(12) 



In either case, if the physical conditions of galaxies at 
higher redshifts are typically well described by compact 
starbursts, the above mentioned energy- loss processes 
will become important relative to synchrotron losses for 
GHz radiating CR electrons. 

To demonstrate this I plot the cooling timescales of CR 
electrons associated with each of the above mentioned 
energy-loss processes versus electron energy in Figure 
[TJ In the left panel of Figure [T] parameters typical for 
normal star-forming galaxies are assumed: B — 10 /iG, 
C/rad = 10~ 12 erg cm -3 , ni SM = 1 cm~ 3 , and ^ osc = 3 kpc. 
Losses due to advection by a galactic wind are ignored. 
At GHz frequencies synchrotron and IC losses are clearly 
the dominant energy-loss terms along with escape. In the 
right panel of Figure [1] the same cooling timescales are 
plotted as in the left panel, however the case of a compact 
starburst having a 1 mG magnetic field is considered. 
Assuming a flux freezing scaling of B oc ^/hism (e.g 



g c 

Ruzmaikin. Sokoloff. fc Shukurov 11198$ iHelou fc Bicav I 



Il993t INiklas fc Beck 1 119971 ICrutcher I Il999f ). it follows 
that nisM = 10 4 cm' 3 . The same ratio of J7s/[/ ra d 
as was used in the left panel is also assumed, setting 
CA-ad = 10~ 8 erg cm~ 3 . Energy losses via advection 
by a starburst-driven galactic wind, having a velocity 
of v w = 300 km s -1 (i.e. typ ical of z ^3 Lyman break 
galaxies; Sha plev et al. 1120031 ) and assuming a disk scale- 
height h — 1 kpc, are also shown. At GHz frequencies 
it is found that synchrotron radiation and IC scatter- 
ing are no longer the dominant energy-loss mechanisms. 
Ionization and bremsstrahlung losses now dominate the 
cooling of CR electrons contributing to the GHz emission 
from galaxies. While such strong fields are present in 
starbursts on small scales, as Zeeman splitting measure- 
ments of OH megamasers in ULIRGs show line-of-sight 
field strengths ranging fro m sa 0.5 — 18 mG in individ- 
ual m asing components (Robishaw, Qua taert. fc Heilesl 
l2008h . given that local starburst galaxies follow the FIR- 
radio correlation, this may suggest that their average 
magnetic fields never reach ~mG strengths since the de- 
pressed synchrotron emission at GHz frequencies should 
yield larger than nominal IR/radio ratios, which is not 
observed. 

2.3. Radio Spectra of Galaxies 

Let us now see how the combination of these energy- 
loss processes translates into the radio spectra of galax- 
ies. For an electron injection spectrum Q(E) — kE~ p , 
the number d ensity of partic les per unit energy is ex- 



partic 

pressed as (e.g lLongair lll994l) 

N(E) 



(P- 1)6(25) 



(13) 
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where b(E) — —dE/dt is the total energy losses for 
electrons including escape. The injected electrons, if 
shock-accelerated by SNRs, will have a spectrum typ- 
ically characterized by Pi n j ~ 2.4. However, as the 
CR electrons interact with the ISM, diffusive losses will 
modify the spectrum, to have an index of p »2.8 (e.g. 
iBiermann fc Strom |[l99l iBecker et al. 112009ft . which is 
assumed for the model. Then, the non-thermal radio 
continuum emission (i.e. the energy losses to CR elec- 
trons by synchrotron radiation relative to the total en- 
ergy losses) is given by 



S„ (xE/r syac N(E), 



(14) 



where emission at frequency v depends on CR electron 
energy E according to Equation [TJ For a constant elec- 
tron temperature, the amount of thermal radio contin- 
uum emission goes as 



oc i/- - 1 . 



(15) 
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Assuming B sw 10 /iG, Ub ~ t4ad> and nisM ~ 1 cm" 
the expected radio continuum spectrum is plotted in 
Figure along with an IR dust emission spectrum 
scaled such that its infrared (IR; 8 — 1000 fim) flux is 
10~ 13 W xnT 2 . The dust emission spect r um is taken 
from the SED libraries of iDale fc Helou I (|2002fl using 
the template which best fit typical star-forming galaxies 
in the Spitzer Infrared Nearby Galaxies Survey (SINGS; 
iKennicutt et al. I [2003) . The thermal radio spectrum 
was then normalized by equating the relations between 
the ioniz ing photon ra te and the thermal radio con- 
tinuu m ([Condon II1992D and IR luminosity (|Kennicutt I 
1998) resulting in 



W ai-' 2 Hz 



VGHz 



6.6 x 10~ 17 
\-o.i 



. w 



( 



(16) 



where T e = 10 4 K is the assumed electron temperature 
for an Hn region. 

Next, the non-thermal component was scaled such that 
the thermal radio fraction is ~10% at 1. 4 GHz, which is 
typically found for star-forming systems (ICondon fc Yin I 
[19901: INiklas. Klein, fc Wielebinlkvl 11997ft . The combi- 
nation of these assumptions naturally leads to the nom- 
inal FIR-radio correlation where 



log 



-Ffir 



3.75 x 10 12 W in- 2 



-log 



4 GHz 



W m^Hz" 1 

, (17) 

is «2.34±0.26 dex ([Yum Reddv. fc Condon |[200l . Re- 
placing the FIR flux with the total IR flux in Equatio n 
[I7]leads to an average q m « 2.64 ±0.26 dex ([Bell Il2003fl . 
For these parameters, it is also found that a least squares 
fit to the non-thermal radio spectrum yields a spectral 
index of «nt ~ 0.8, simila r to what is observed in normal 
star-f orming galaxies (e.g. INiklas. Klein, fc Wielebinskv"! 
[1991 . 

To keep a fixed ratio between the FIR and radio emis- 
sion of galaxies requires that the ratio of synchrotron to 
total energy losses for CR electrons to be nearly con- 
stant among each system. Or, in the case of normal 
star- forming galaxies, where CR electron energy losses 
are dominated by synchrotron losses and IC scattering 



2.5 - 



2.0 - 













0.5 GHz 1 




\ 


0.3 GHz : 






'. . £3 — 1 UU jui\J 






1 ■ ■ - 50 fiG 




B=10 fiG 





Fig. 4. — The ratio of observed-frame 0.3 and 0.5 GHz to 1.4 GHz 
flux densities plotted as a function of redshift in blue and red, re- 
spectively. These estimates include the expected suppression of 
non-thermal radio continuum as IC scattering of th e CMB become 
increasingly important at higher rcdshifts (sec i|3,ll for details about 
the calculations.) Different cases for the intrinsic magnetic field 
strength of galaxies are considered: 10 /iG (solid lines), 50 /iG 
(dashed lines), and 100 ^G (dottted lines). Unless magnetic field 
strengths are very strong (> 10 fiG), only marginal sensitivity gains 
exist by making lower frequency radio observations to detect galax- 
ies at z > 2. 

off a galaxy's interstellar radiation field (ISRF), the FIR- 
radio correlation is preserved under the condition where 



U B 



ad 



>i, 



(18) 



where C/cmb is the radiation field energy density of the 
cosmic microwave background (CMB). In the local Uni- 
verse, IC losses off of the CMB are negligible to those 
occurring off the ISRF within a galaxy, however this is 
not necessarily the case at higher redshifts. 

3. RESULTS: EXPECTATIONS FOR DEEP RADIO 
CONTINUUM SURVEYS 

I now explore implications of the above physical ar- 
guments for evolution in the observed radio continuum 
emission from galaxies at increasing redshifts. I also re- 
port on how this relates to an expected evolution in the 
FIR-radio correlation; at present, it is assumed that any 
change in the FIR-radio correlation arises solely from 
changes to the radio continuum emission of galaxies. 
Variations arising from differences in the IR emission 
of galaxies are discussed below in ^4.4. 21 Implications 
for how various galaxy properties, including their star 
formation activity and magnetic field properties, can be 
inferred from the discrepancies between expected and ob- 
served IR/radio ratios are discussed in §4. 

3.1. Suppression of Non- Thermal Emission by the 

CMB 

As already shown, the FIR-radio correlation relies on 
a fixed ratio between synchrotron and the total energy 
losses of CR electrons. At z = 0, Uqmb ~ 4.2 x 
10~ 13 erg cm~ 3 which is significantly smaller than the 
radiation field energy density of the Milky Way (i.e. 
Umw ~ 10~ 12 erg cm" 3 ). Thus, CR electron energy 
losses from IC scattering off the CMB are negligible at 
low redshifts. However, Uqmb oc (1 + z) 4 making such 
losses increasingly important with redshift. For instance, 
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by z ~ 3, E/cmb ~ 1-1 x 10 -10 erg cm -3 ; equating this 
to Ub = B 2 /(8ir) results in a corresponding magnetic 
field strength of ~50 /iG, nearly an order of magnitude 
larger than the ambient field strength in the Solar Neigh- 
borhood. Consequently the non-thermal component of 
a galaxy's radio continuum emission will be increasingly 
suppressed with increasing redshift, eventually resulting 
in only the thermal component being detectable. 

Assuming an intri nsic FIR-ra dio correlation for galax- 
ies with <7ir ~ 2.64 (jBell I [2003h . the expected observed- 
frame 1.4 GHz flux density for star- forming galaxies hav- 
ing a range of IR (8 — 1000 /im) luminosities as a function 
of redshift are estimated (Figure[3]). These estimates rely 
on the results of the previous section where, on a theoreti- 
cal basis, a realistic radio spectrum was constructed. It is 
assumed that the radio continuum emission is comprised 
of two components, thermal (free-free) and non-thermal 
(synchrotron) emission, both of which can be expressed 
as power-laws such that S u cx v~ a , where «t ~ 0.1 for 
the thermal (free-free) component and q;nt ~ 0.8 for the 
non-thermal component (Table [TJ). A thermal fraction of 
~10% at 1.4 GHz is also assumed. And a (1 + z)^ 1 cor- 
rection factor to account for bandwidth compression as a 
function of redshift is included. Increased energy losses 
to CR electrons from IC scattering as C/cmb increases 
with redshift are explicitly included. Figure [3] illustrates 
two cases per IR luminosity in which the galaxy has a 
10 and 100 /iG field, indicated by the solid and dotted 
lines, respectively. The case of a moderately bright LIRG 
(Lir ~ 3 x 10 11 Lq) is highlighted in red, and the case 
of a 50 /iG field strength shown using a dashed line. 

For the case of a moderately bright LIRG, it is clearly 
shown that its detection in the radio continuum may sig- 
nificantly rely on the strength of its magnetic field be- 
tween 1 < z < 6 due to the suppression of its synchrotron 
emission by increased IC losses off of the CMB. For ex- 
ample, a galaxy having B = 10 /iG at z ~ 3 will also 
be detected at z ~ 4.5 if its magnetic field were 100 /iG. 
At higher redshifts it appears that the strength of the 
magnetic field plays less of a role since IC losses off of 
the CMB begin to completely suppress a galaxy's non- 
thermal emission, making only its thermal radio emission 
detectable. It is also shown that a moderately bright 
LIRG at z ~ 10 should have a 1.4 GHz flux density of 
~40 njy almost entirely arising from free-fr ee emission. 
Using the star formation rate conversion of iKennicutt I 
(1998), this implies that all galaxies forming stars at 
a rate of > 50 M Q yr _1 should be detectable out to 
z ~ 10 if such a population of galaxies exists at these 
early epochs and this sensitivity is achievable. 

Due to the relatively IR/radio negative spectral slope 
of a galaxy's non-thermal radio continuum emission, the 
case of observing at lower frequencies is considered to see 
what sensitivity gains might exist for detecting galaxies 
at high-z. The expected observed- frame flux densities at 
0.3 and 0.5 GHz are estimated. The ratio between the 
observed-frame 0.3 and 0.5 to 1.4 GHz flux densities as a 
function of redshift is plotted in Figure [4j The cases for 
different intrinsic B- field strengths among galaxies (i.e. 
10, 50, and 100 /iG) are again considered. For galaxies 
having magnetic field strengths of ~10 /iG it appears that 
there is little advantage for surveys at lower frequencies 
since, by z ~ 2.5, sensitivity gains are less than ~50%. 



If a galaxy's magnetic field is quite strong, however, be- 
ing ^100 /iG, observations at 0.3 GHz should be more 
sensitive to normal star-forming galaxies than 1.4 GHz 
observations by a factor of >2 out to z ~ 6. While 
the possibility of increased magnetic field strengths in 
galaxies at higher redshifts is considered, there is not a 
strong physical basis to expect this as the leading the- 
ory of magnetic field generation, the mean-field dynamo 
model, predicts large-scale regular magnetic fields to be 
weak er in the past rather than stronger (e.g. lBeck et al. I 
Il99l . 

3.2. Expected Evolution in the FIR-Radio Correlation 

with Redshift from CMB Effects 

As already shown, the non-thermal component of a 
galaxy's radio continuum emission is increasingly sup- 
pressed as a function of redshift due to the increasing 
importance of IC scattering off of the CMB. This ef- 
fect should be directly reflected in changes to a galaxy's 
IR/radio ratio. Assuming that the intrinsic IR/radio ra- 
tio of star-forming galaxies is not unique to z ~ sys- 
tems, the FIR-radio correlation should evolve in a very 
predictable way as a function of redshift arising from the 
increase in the energy density of the CMB. This is illus- 
trated in Figure [5] When the non-thermal emission is 
completely suppressed by IC scattering off of the CMB, 
observed IR/radio ratios are expected to simply approach 
the ratio between the FIR and thermal radio continuum. 

As already stated, different CR electron cooling pro- 
cesses may become increasingly more important than 
synchrotron radiation in galaxies hosting strong star- 
bursts. For such cases, the non-thermal component 
may diminish even more quickly than in the calculation 
provided here, resulting in observed IR/radio ratios for 
which only the thermal radio continuum emission is de- 
tected (i.e. qin ratios which are a factor of ~10 higher 
than the nominal value). However, it is also possible 
that other physical processes could work in the opposite 
direction to increase the synchrotron emissivity, such as 
an increase to the acceleration efficiency of CR electrons 
within SNRs, as well as additional synchrotron emission 
arising from an increase to secondary positrons and elec- 
trons. These possibilities, though speculative since such 
local starbursts do not seem to exhibit such variations in 
their IR/radio ratios, are discussed in 8 )4.4.31 and 14.4.41 
respectively. 

3.3. Scientific Goals and Technical Specifications for 
Deep Radio Continuum Surveys using the SKA 

Detecting moderately bright star-forming galaxies at 
redshifts beyond z ~ 3 requires extremely sensitive con- 
tinuum imaging only expected with the advent of the 
SKA. Figure [3] shows the achieved point-source sensi- 
tivities for a number of existing 1.4 GHz radio con- 
tinuu m surveys including: VLA-COSMOS ( S chinnerer I 
[20071 ctrms ~ 50 fjjy), GOODS-N (G. Morrison 2009, in 
prepar ation, 5ctrmk ~ 20 /xJy and the Deep SWIRE 
Field (|Owen fc Morrison I 120081 5cr RM s «13.5 /iJy). 
The expected point-source sensitivity of the EVLA af- 
ter 300 hr (~1 Ms) is also plotted (5ctrms ~ 1-4 /zJy). 
As shown, even the EVLA will not be able to de- 
tect individual M 82-like objects, having IR luminosities 
-4 x 10 10 L Q , at z > 1. 
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Fig. 5. — The expected rest-frame IR/radio ratios for a galaxy 
having qui = 2.64 at z = as IC losses off of the CMB become 
increasingly important as a function of redshift. Each track corre- 
sponds to a different internal magnetic field strength for the galaxy. 
As Ucmb 3> Ub, the IR/radio ratio approaches the limit where 
only thermal (free-free) radiation contributes to the observed ra- 
dio continuum emission. The average local qi^ values (2.64 dex; 
dot-dashed line) and the ±1 a scatter ( dotted-Yme) are shown. 

As discussed in ij3.ll a 5c point-source sensitivity of 
~40 nJy is required for the detection of a moderately 
bright LIRG (i.e. Lir ss 3 x 10 11 £©), forming stars at 
~50 Mq yr _1 at z ~ 10. However, it is unclear as to 
whether there will be a significant population of dusty 
star-forming galaxies at these redshifts. At present, only 
rest- frame UV studies relying on the "dropout" tech- 
nique are able to probe normal galaxies at such high 
redshifts. In Figure [3] the estimated observed-frame 
1.4 GHz flux densities, based on the results from the UV 
(w 1600 A) l uminosity function work of iBouwens et alH 
(20071 l2008h . are shown. These authors use space- and 
ground-based data over both Hubble ultra-deep fields and 
GOODS fields to search for star-forming galaxies at z ~4, 
5, 6, 7, and 9. The UV luminosities, ov er the full range 
per re dshift bin as gi ven in Table 5 o f IB ouwens et al. I 
i|2007h and Table 3 of IBouwens et al. I (12008}) . wer e cor- 
rected for dust following IBouwens et al. I (f2007) such 
that v4i 6 oo = 1.1 mag, 0.6 mag, and 0.4 mag for z ~ 4, 
z ~ 5, and z > 6 galaxies, respectively. Corresponding 
observed-frame radio flux densities were then estimated 
by first converting the extinction -corrected UV lum inosi- 
ties to IR luminosities following IKennicutt I (|1998f ). and 
then using the FIR-radio correlation assuming a galaxy 
magnetic field strength of 10 fiG. 

These considerations suggest that the criterion of 
a point-source sensitivity of «40 nJy will only skim 
the brightest sources among this population of galax- 
ies. Beyond z > 7, detection of these galaxies will be 
marginal at best. However, it is possible that the as- 
sumed dust corrections significantly underestimate the 
true star formation rates of these sources as they are 
based on an empirical relation between the UV contin- 
uum slope and extinction in local starbur sting galaxies 
(e.g. IMeurer. Heckman. & Calzetti I [19991) . While this 
prescriptio n works well for z ~ 2 LIRGs older than 
~100 Myr (|Reddv et al. II2006D . it has al so been shown to 
fail for many other high-z systems (e.g. iGoldader et al. I 
l200l IReddv et al. 1120061: IMurphv et al. I I2009al) . Ym- 
thermore, the 4 < z < 9 UV luminosity function work 



suffers dramatically from the effects of cosmic variance. 
In any case, it appears quite advantageous for the SKA 
sensitivity requirement be a factor of ~2 better, namely 
a 5a point-source sensitivity of « 20 nJy (i.e. crms ~ 
4nJy). 

The sensitivity for the SKA is approximated us- 
ing the radiometer equation such that ctrms ~ 
2k B /y/2 BW ^nt (Aeff/Tsys)- 1 where k B is the Boltz- 
mann constant, BW is the bandwidth, i int is the inte- 
gration time, A c s is the effective collecting area of the 
array, and T sys is the system temperature. Assuming 
a bandwidth of ~1 GHz (comparable to what existing 
telescopes are delivering) and a reasonably long integra- 
tion (i.e. 300 hr), I find a requirement for the SKA of 
A c g/T sys ~ 15000 m 2 K _1 to reach this sensitivity such 
that, 

(<trms\ „ 4 /BW\ - 1/2 f t iat \ - 1/2 f A cff /T sys \ - 1 

\nJy)~ \GRzJ V 300 W V 15000 m2 K_1 / 

(19) 

With these specifications, the SKA will be sensitive to all 
galaxies forming stars at a rate of > 25 Mq yr _1 . The 
5ctrms ~ 20 nJy is also plotted in Figure [3l 

3.3.1. Technical Challenges for Deep SKA Continuum 

Fields 

In order to achieve such highly sensitive radio imagery, 
two major technical challenges must be considered in 
the SKA design to ensure that systematic noise does 
not limit the usability of such imaging; natural confu- 
sion from overlapping sources in the field and insufficient 
dynamic range to identify s ources properly. Both effects 
were recently considered bv lCondoiil (|2009l ). who found 
that they should not limit deep continuum surveys using 
the SKA at nJy depths. 

While the angular size distribution of the sub-^Jy 
population remains uncertain, most faint radio sources 
are likely coextensive with faint star-forming galaxies at 
z 1. These sources are therefore probably smaller than 
~l"setting the natural confusion limit to be the flux den- 
sity corresponding to ^30 arcsec~ 2 per source. Accord- 
ingly, the synthesized beam of the SKA should be 6$ < 1" 
to keep instrumental confusion below the natural confu- 
sion limit , requ ir ing baselines up to ~50. km at 1.4 GHz 
(|Condon 1 120091 ). iWindhorst elaTl ([200l) also estimate 
the expected 1.4 GHz source sizes at 10 — 100 nJy by 
taking the galaxy disk size distributions measured from 
deep HST imaging, where the object density is in excess 
of 10 6 objects deg -2 . Their findings suggest that the 
synthesized beam of the SKA be even sharper, reaching 
#S ^0'.'3, to avoid instrumental confusion and requiring 
baselines up to ~150 km at 1.4 GHz. The 1.4 GHz natu- 
ral confusion is es timated to be a c < 3 nJy, though with 
large uncertainty ICondonl (|2009l ). suggesting that from 
the point-source sensitivity calculation above, the SKA 
should not be limited by natural confusion. 

Coupled with the requirement of exquisite sensitivity 
is achieving an unprecedented dynamic range to actually 
make use of such deep imagery. Typically, the dynamic 
range is defined as D = S max / 'ctrmS: the ratio of the 
strongest peak (5 max ) to RMS fluctuations in regions not 
containing sources, and describes an image dominated by 
an extremely strong point source. At present, the highest 
dynamic range achieved has been ^60 dB. For example, a 
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TABLE 1 

Assumptions for Expected Radio 
Flux Densities 



Parameter 


Value 


<?IR 


2.64 




0.1 


a NT t 


0.8 


f lA GHz 

It 


0.10 


Ucmb(z = 0) 


4.2 X 10 _13 ergs cm -3 



Note. — t; The thermal and non- 
thermal radio continuum is assumed to 
follow power-laws of the form Si, <x 
v~ a , where «t and aNT are the as- 
sumed thermal and non-thermal spec- 
tral indices. 



dynamic range of ~60 dB was obtained with Westerbork 
Synthesis Radio Telescope (WSRT) observations over a 
10' region containing a single bright source (i.e. 3C 147). 
At the desired nJy sensitivity, the large primary beam 
of the SKA will contain many sources having strengths 
larger than that of the radiometer noise. Thus, the dy- 
namic range required by the SKA can be estimated using 
the typical flux density for the strongest source in a field, 
which for a primary bea m f2b ~ 1 deg 2 , is ^0.16 Jy at 
1.4 GHz (Condon 2009). Using the required RMS noise 
of orms ~ 4 nJy, this leads to a dynamic range require- 
ment of D > 4 x 10 7 = 76 dB. Reaching this dynamic 
range will likely require, at minimum, antennas having 
extremely stable beams and pointing capabilities. Such 
high precision antennas should function well at frequen- 
cies >10 GHz, which is the current frequency cutoff of 
the SKA Phase- 1 and Phase-2. 

Table [5] summarizes the requirements on the SKA in 
order that it be able to detect galaxies forming stars with 
rates >25 M© yr _1 at all redshifts where such galax- 
ies may be expected. As such, the SKA would there- 
fore be sensitive to a significant number of the high-z 
g alaxies reported in the UV lu minosity function studies 
of lBouwens et al. I (|2007ll2008f) . However, there are still 
large uncertainties in the current understanding of the 
sub-/ijy population of radio sources, and achieving such a 
high imaging dynamic range represents a significant tech- 
nical challenge. The current generation of SKA pathfind- 
ers (e.g. EVLA, ASKAP, MeerKAT, and APERTIF) 
could usefully reduce the uncertainties in both the sub- 
/zJy radio source population(s) and high dynamic range 
imaging. 

Finally, as the observed radio continuum emission from 
normal star-forming galaxies at high redshift will likely 
be dominated by their thermal emission, this analysis 
does not require that dust be present in these galaxies 
for these sensitivity calculations to be valid. Specifically, 
this analysis requires their bolometric luminosities, and 
not necessarily their IR luminosities, be > 10 11 L© for 
detection. Thus, radio continuum emission from galax- 
ies becomes the only way to find and quantify the star 
formation properties of galaxies unbiased by dust. The 
importance of this fact for studies of the cosmic star for- 
mation history are discussed in i j4.31 

4. DISCUSSION 



TABLE 2 

SKA Deep Continuum Survey Assumptions 



Parameter 


Value 


Target Star Formation Rate 


25 M Q yr" 1 


Bandwidth (BW) 


1 GHz 


Angular Resolution (8$) 


< tf'3 - 1" 


Primary Beam (fit,) 


"» 1 deg 2 


Imaging Dynamic Range (D) 


> 76 dB 


Integration Time (ti nt ) 


300 hr (~1 Ms) 


SKA Sensitivity (A eS /T sys ) 


15000 m z K _i 


Image RMS Noise Level (optMs) 


«4 nJy 



I will now discuss how various galaxy properties can 
be inferred from discrepancies between expected and ob- 
served IR/radio ratios. I will also discuss what this 
means for the prospect of using deep radio continuum 
surveys, both at high and low frequencies, to study star 
formation over cosmic time using the SKA. 

4.1. Searching for the Emergence of Large- Scale 
Magnetic Fields 

One of the primary energetic constituents of the ISM is 
its magnetic field. Its energy density is equal to that of a 
galaxy's radiation field, CRs, and turbulent gas motions, 
yet its origin and role in shaping galaxy evolution, in par- 
ticular through the possible r egulation or star formation , 
remains highly uncertain (e.g. lFerriere 1120011 ICox 1 12005). 
At present the leading theory for large-scale magnetic 
field generation i s the mean-field dynamo model (e.g. 
iBeck et al. 1 119961 ). In this paradigm galaxy-scale mag- 
netic fields are built up over time, suggesting that field 
strengths should be typically much weaker in the past. 
The exact details on how the amplification of weak, seed 
magnetic fields in the halos of protogalaxies eventually 
translates into galaxy-scale fields is uncertain. Does the 
first epoch of star formation quickly initiate a dynamo 
in a galaxy? Does the field form during the in-fall of gas 
prior to star formation? Recent theoretical work suggests 
that weak, seed magnetic fields may be quickly amplified 
by the turbulent (small-scale) dynam o, yielding strong, 
(~20 /zG) turbulent fields by z ~ 10 {Arshakian et al. 
2009). These fields then serve as a seed for the mean- 
field (large-scale) dynamo, resulting in regular fields of 
/iG strengths having a few kpc coherence within ^2 Gyr 
(i.e. by z ~ 3). Such a scenario is consistent with indirect 
evidence from Faraday rotation studies of distant radio 
sources which clai m /iG fields to be in place for normal 
galax ies at z ~ 1.5 (Kr onberg et al. Il2008l ; iBernet et al. I 
120081) . 

Since the non-thermal radio emission from galaxies, 
and therefore the linearity of the FIR-radio correlation, 
depends critically on the magnetic field strength of a 
galaxy, deviations from nominal ratios as a function of 
redshift may help to illuminate the presence and char- 
acteristics of magnetic fields in galaxies at early epochs. 
However, given that regular magnetic fields are thought 
to weaker in the past, as suggested by the mean- field dy- 
namo model, maintaining a constant IR/radio ratio with 
increasing redshift by increasing magnetic field strength 
with redshift is somewhat contrived. If the FIR-radio 
correlation is found to hold for normal galaxies at z >3, 
either the the turbulent dynamo and associated small- 
scale magnetic fields play an important role at early 
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epochs, or our physical picture for the correlation needs 
revision. This assumes of course that the existence of an 
AGN in these galaxies, whose radio emission may for- 
tuitously bring the IR/radio ratio to near or below the 
local correlation value, can be ruled out. Such situations 
are discussed below in §32] where I compare observations 
with these expectations. 

In the absence of strong magnetic fields at high-z, the 
IR/radio ratios of galaxies should be large, similar to 
what is observed for 'nascent' starburs ts. These sys- 
tems are rare (<1%) in the local Universe (Roussel et al. I 
2003), and are thought to be experiencing the onset of 
an episode of star formation after a long (wlO 8 yr) pe- 
riod of quiescence. Since OB stars will always ionize the 
nearby gas, generating free-free emission, these galaxies 
may only be identifiable in the radio by their thermal 
radio continuum emission. 

4.2. Current Observational Results 

Using existing facilities a number of deep FIR, 
submm, and mm surveys have detected galaxies at 
fairly high rcdshfits. Sources detected in the FIR 
and radio from a number of these studies are in- 
cluded in Figure [6] with the highest re dshift sources 
being the z ~ 6 quasar s detec ted by iCarilli et al. I 
(2004) and iBeelen et al~l (|2006f ). In cases where 
only a FIR measurement is provided, the assump- 
tion that IR(8 - 1000 /j,m) « 2 x FIR(42 - 122 fim) is 
made, which is exactly the difference between the 
average FIR/radio and IR/ r adio ratio s repo rted by 
lYun. Reddv. & Condon I (pOOl andEeh] ([2003h . respec- 
tively. The average and ± ler lines for the local IR- radio 
corre lation are drawn (i.e. <7ir ~ 2.64 ± 0.26 dex; iBell I 
20.03). As in Figure[SJ the expected locations of a typical 
star forming galaxy as IC losses of the CMB become in- 
creasingly important relative to the internal synchrotron 
losses of CR electrons for various magnetic field strengths 
are overplotted. Since so few sources are detected at 
z > 3, a logarithmic scaling is used for the rcdhift axis 
to better identify any trends. 

4.2.1. Submillimeter Galaxies 

Most of the galaxies out to z ~ 1 (e.g. the 
70 /im detected sources in; iFraver et al. I [20061: 
IE. J. Murphy et al. 112009b. in preparation! ) tend to have 
IR/radio ratios similar to the local value, which is consis- 
tent with the expectations presented here since the mag- 
netic fields in such galaxies need not be excessively large 
to compete with IC losses off of the CMB. Once beyond 
z *~ 1 there is an increasingly larger number of galaxies 
having IR/radio ratios below the canonical value, rather 
than larger per the suggested expectations. The major- 
ity of these sources are submillimeter galaxies (SMGs). 
While typically detected in the hard-band (2 — 8 keV) X- 
rays, thus indicating the presence of an AGN, this result 
is somewhat surprising since SMGs are thought to be pri- 
marily powered by star form ation (e.g. lAlexander et al. I 
l2005t Murphy et al. Il2009al). 

The sampl es of IKovacs et al. I (|2006| ) and 
iValiante et al. I (|2007f ) consist entirely of SMGs 
and the filled t riangle s indicate the SMGs included in 
IMurphv et aT~l (|2009a). The mean redshi ft among these 
three samples of SMGs is z ss 2.17 ±0.65. IKovacs et al. I 



( 2006) report IR/radio ratios for their sa mple of SMGs to 
be syst ematically below the local value. IMurphv et al. I 
(2009a) also found that the SMGs in their sample to 
have IR/radio ratios which are, on average, a factor 
of ^3.5 lower than the canoni cal ratio. The mean 
IR/radio ratio for the sample of IValiante et al.~l (2007) 
is ~ 2.27 ± 0.16 dex, a factor of ~2 below the local value 
of 2.64 dex. The average IR/radio ratio among all 38 of 
these z ~2.2 SMGs is q m = 2.26 ± 0.41 dex, a factor of 
~2.4 times smaller that the local value. A number of 
possible explanations for why these sys tems depart for 
the FIR -radio correlation are included in IMurphv et al. 1 
( 2009a) , but a definitive conclusion has not been drawn. 

Recently, a number of higher redshift SMGs 
have been 
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identifi ed dCapak et al. I 
2009allbl: iCoppin et al. ' 



2008 



2009 



V. Smolcic et al. 112009. in preparation!: I A. Karim et al. I 
2009. in preparation! : iKnudsen et al. II200^ K The aver- 
age redshift among these 8 objects is z ss 4.40 ± 0.33. 
Similar to the z ~ 2.2 SMGs, the mean IR/radio ratio 
of these z ~ 4.4 SMGs is q m = 2.20 ± 0.27 dex, a factor 
of ~2.8 times smaller that the local value. This result is 
somewhat surprising since the expected increase in the 
IR/radio ratios between a z ~ 2.2 and z ~ 4.4 galaxy 
due to the increase in IC losses to CR electrons from the 
CMB is a fact or of ~3 for a ~35 /iG field, the average 
value for the IMurphv et al~l (|2009af) SMGs using the 
revise d minimum energy formula of iBeck &: Krausel 

While it has been suggested that the properties of the 
z ~ 4.4 SMGs are similar to the z ~ 2.2 SMGs (e.g. 
iDaddi et al. Il2009al) . the consistency in their IR/radio 
ratios argues the opposite, suggesting that their star for- 
mation driven IR and/or radio continuum emitting prop- 
erties are sufficiently different to account for the large dis- 
crepancy arising from extra cooling of CR electrons by 
IC scattering off the CMB. Alternatively, this result may 
also simply suggest that, unlike their FIR output, the 
radio emission from both populations of SMGs is in fact 
dominated by AGN, or that their magnetic fields are sig- 
nificantly strong (i.e. > 300 fxG) such that Ub ^> Ucmb 
at these redshifts. The latter explanation is tough to 
reconcile since, in the presence of such strong fields, CR 
electrons emitting at GHz frequencies will be cooled pri- 
marily by ionization and bremsstrahlung losses rather 
than by synchrotron radiation (see §2.2), unless some 
other physics is at work. At present it is difficult to 
make any strong conclusions from this comparison since 
it suffers from low number statistics; very few SMGs that 
have been identified at z > 4. 

4.2.2. Radio Quiet QSOs at High-z 

The samples of lCarilli et al. I (|2004h and IBeelen et al. I 
(2006) consist entirely of radio-quiet quasars QSOs, a 
total of 6. The highest redshift QSO, J114816+525150 
(hereafter, J1148+5251), is identified in Figure[B]and ap- 
pears to have an IR/radio ratio which is amazingly con- 
sistent with the local value. In fact, 4/6 (i.e. ~66%) of 
the QSOs have IR/radio ratios that are within la of the 
canonical ratio. While these observations are suggestive 
that the energetics of these objects is dominated by star 
formation, such a conclusion is not so clear. 

These galaxies may lie near the nominal IR/radio ratio 
due to excess emission from their AGN. Even if all of the 
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Fig. 6. — The same as Figure [5] except that I have plotted actual data from a number of studies in the literature. A logarith- 
mic re^sliift_scale is_used since there are so few points at z > 3. F i lled in symbo l s indic ate that an object is an SMG, i.e. the en- 
tire IKovacs et al. I 112003). I Valiante et al I ll2007l)ICapak et al nl200g|). iDaddi et al I H2009al lH), ICoppin et al~l (I2009T) . IV. Smolcic et al. I 
igOOg jin prep aration). A. Kari met al. I 112009, in preparation !) . and Knudsen et al. (2009j) samples, as well as half of the IMurphv et al. I 
Il2*0*0*9*a!r sample. The lower limit (uvvward-arrow) for the iKnudsen et aT!Tll2*0*0*9l) data point is due to only having an upper limit at 
1.4 GHz. The radio-quiet quasar J114816+525150 is labeled; at a redshift of z = 6.64 it lies at ~1/16 of the Universe's current age. Rather 
than seeing galaxies tend toward higher IR/radio ratios per expectations given that IC losses of the CMB become increasingly important 
with increasing redshift, many galaxies tend to have IR/radio ratios lower than the local value (i.e the SMGs). 



synchrotron emission associated with ongoing star for- 
mation in these galaxies is suppressed due to other CR 
electron cooling processes competing successfully with re- 
spect to synchrotron radiation, the excess emission from 
an AGN may be enough to fortuitously put a galaxy 
back on the local correlation. This scenario is likely the 
simplest way to explain how QSO's at these redshifts 
have IR/radio ratios consistent with normal star-forming 
galaxies in the local Universe. 

Alternatively, these objects would require magnetic 
field strengths nearing >500 fj,G levels to remain on 
the local FIR-radio correlation if powered by star for- 
mation alone. Taking Jl 148+5251 as an example; the 
galaxy is being observed <1 Gyr after the Big Bang. 
While it has been suggested that strong (^20 //G) 
turbulent fields may be present by z ~ 10, it is 
unclear how these seed fields could be amplified so 
rapidly and form galaxy-wide large-scale regular fields 
of > 200 — 500 /iG. If global star formation initiates a 
large-scale dynamo which acts to regulate the magnetic 
fields of galaxies fe.g.lMichel fe Yahilll 19731 : [Reese II1987K 
iRuzmaikin. Sokoloff. fc Shukurov 1 119881 ). the timescale 
for this must occur rapidly as star formation is only 
~400 Myr old assuming that star formation began 



around z ~ 10. The requirement for such a rapid am- 
plification phase suggests that the large-scale field may 
begin to form before the initiation of star formation, dur- 
ing the in-fall of gas and dust (e.g. [Pudritz fe Silk 1 1 1 989l: 
Beck et al. 1 119941: iHelou fc Bicav Ml 9931: iKulsrud et al. I 
19971) . However, these strong magnetic fields would then 
have to decay considerably over cosmic time to recover 
the typical strengths observed for galaxies in the local 
Universe. Detections of less exotic systems at these red- 
shifts in the FIR and radio (e.g. using ALMA and the 
SKA) will be necessary to reconcile these outstanding 
questions. In fact, a comparison of the IR/radio ratios 
for these high-z galaxies will likely be one of the only 
methods to test for the existence and strength of mag- 
netic fields in the early Universe while trying to assess 
what role they may play in star formation. 

4.2.3. Lyman Break Galaxies at z ~ 3 

ICarilli et al. I (|2008l ) recently analyzed the radio prop- 
erties for a large sample of Lyman break galaxies (LBGs) 
at z ~ 3. Using a stacking analysis, and assuming 
the local FIR-radio correlation, they find that the ra- 
dio (1.4 GHz) derived star formation rate is a factor of 
1.8 ± 0.4 times larger than that from observed UV mea- 
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surements. Assuming that this difference arises from at- 
tenuation by dust at UV wavelengths, this value is signif- 
icantly less than the factor of ~5 normally assumed for 
LBGs. They suggest that this discrepancy may be the re- 
sult of depressed radio emission resulting from increased 
IC cooling of CR electrons off of the CMB, consistent 
with the analysis presented here. 

The stacked 1.4 GHz flux density reported by 
iCarilli et al~l (|2008f ) is 0.90±0.25 /xJy, which corresponds 
to a rest-frame 1.4 GHz flux de nsity of 0.68 /iJ y assum- 
ing a radio spectral index of 0.8 (|Condon Ill992f). Taking, 
the t ypical size of an LBG to be <1" (e.g. iGiavalisco I 
21)021 results in a minimum e nergy magnetic field of 
> 8 fiG (|Beck fc Krause I [20051 . The minimum energy 
magnetic field is relatively insensitive to the radio surface 

brightness (i.e. i? m in oc ll^ 3+aS> ); even if the reported 
1.4 GHz flux density were under-predicted by a factor 
of 10, the corresponding magnetic field strength would 
only increase by a factor of ssl.8. The increased energy 
loss to CR electrons from IC scattering off of the CMB 
at z ~ 3, assuming equipartition between the galaxy's 
magnetic and radiation field energy densities, will de- 
crease the observed 1.4 GHz flux density by a factor of 
~3.5 which is easily able to account for the discrepancy 
between dust attenuation factor suggeste d by the radio 
(1.4 G Hz) derived star formation rate of ICarilli et al. I 
(2008) and the factor of ^5 normally assumed for LBGs. 
If the intrinsic magnetic fields of LBGs are as large as the 
minimum energy fields for SMGs (i.e. ~ 35 fiG,) then the 
observed non-thermal radio continuum emission should 
still be suppressed by a factor of ~2. 

4.3. Implications for Radio-Based Studies of Star 
Formation at High Redshifts 

While I have looked to see how the GHz emission from 
galaxies will be affected by increased CR electron cooling 
as the result of IC scattering off of the CMB at high red- 
shifts, I have yet to look at what this may do for higher 
frequency radio observations. I now consider the poten- 
tial usefulness of low (<2 GHz) and high (>10 GHz) fre- 
quency radio observations for measuring star formation 
over cosmic time. 

4.3.1. At Low Frequencies: <2 GHz 

As discussed above, studies of star formation at low 
frequencies using the SKA will be complicated by the 
expected evolution in the FIR-radio correlation with in- 
creasing redshift. When coupled with the additional 
emission from AGN, which may appear to place galax- 
ies back on the FIR-radio correlation, it appears that 
using GHz observations to estimate the cosmic star for- 
mation history may lead to increasingly unreliable esti- 
mates of star formation rates with increasing redshift. 
Galaxies which host compact starbursts may have large 
magnetic fields and gas densities in which case electrons 
contributing to GHz emission will lose their energy pri- 
marily to ionization and bremsstrahlung processes rather 
than synchrotron radiation. These galaxies may move off 
the correlation by some fixed amount depending on how 
the various energy-loss terms compete, resulting in incor- 
rect radio-based estimates of star formation rates using 
locally-derived calibrations. In conclusion, it seems that 
interpreting the non-thermal GHz measurements of high- 
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Fig. 7. — The same as Figure [3] except that the expected emis- 
sion at 10 GHz is considered. The expected sensitivity of the 
SKA at 10 GHz is also assumed identical to the 1.4 GHz sensi- 
tivity. Since the non-thermal contribution is almost completely 
suppressed at high redshift, the 10 and 1.4 GHz emission should 
be detectable from the same population of galaxies at these epochs. 
The expected flux density range for galaxies incl uded in the z ~ 
4, 5, 6, 7, an d 9 UV luminosity function studies of Bouwcns et al. 
1)20071 . |20Q8) is again given by the shaded region. 
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Fig. 8. — The expected thermal fraction for observed-frame 1.4 
(blue) and 10 GHz (red) observations of normal star-forming galax- 
ies vs. redshift. Intrinsic magnetic fields of 10 (solid-lines), 50 
(dashed-lines) , and 100 fiG (dotted-lines) are shown. By observing 
at 10 GHz, the observed radio emission is dominated by thermal 
processes, making it an ideal measure for the massive star forma- 
tion rate of galaxies, beyond a redshift of z ~ 2 even for magnetic 
field strengths of ~100 fiG. 

z galaxies properly may prove difficult. 

4.3.2. At High Frequencies: >10 GHz 

While deep continuum studies aimed at measuring the 
star formation properties of high-z galaxies at low fre- 
quencies (i.e. <2 GHz) are not straightforward due 
to large variations in the amount of non-thermal radio 
emission, such studies at slightly higher frequencies (i.e. 
>10 GHz) are likely to be much more promising. In 
fact, it can be argued that deep radio continuum stud- 
ies at these frequencies will provide the most accurate 
estimates for the star formation rate of high-z galaxies. 
So long as average galaxy magnetic fields do not exceed 
mG strengths, which seems unlikely given that at such 
high redshifts it is unclear how such strong fields could be 
seeded, the observed radio emission at >10 GHz should 
be dominated by thermal (free-free) radiation directly 
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associated with Lyman continuum photons emitted by 
newly formed massive stars. 

In Figure [7] the expected observed- frame 10 GHz con- 
tinuum flux densities of galaxies are plotted for a range 
of IR luminosities under the same assumptions that went 
into creating Figure [3] (see Table H}. Overplotted is the 
expected range in observ ed-frame 10 GHz flux densit y 
for the sources studied by iBouwens et al. I (|2007l I2008T ). 
It is shown that even at 10 GHz the same population 
of LIRGs, if they exist, should be detectable at all red- 
shifts. Furthermore, it appears that the internal mag- 
netic strength of a galaxy does not play as signifiant a 
role in a source's detectability since the thermal fraction 
at these frequencies is so large, and becomes increasingly 
so with redshift. This is illustrated in Figure [8] where the 
expected thermal fraction for observed 1.4 and 10 GHz 
flux densities is plotted as a function of redshift for galax- 
ies having 10, 50 and 100 /iG fields. 

Observations at 10 GHz are dominated by thermal 
emission for B < 100 /iG beyond a redshift of ~2. The 
same cannot be said for 1.4 GHz observations until the 
galaxy population being observed lies at a redshift of 
z >5. At 10 GHz, the radio emission from z > 5 galax- 
ies is likely to be >80% thermal emission. Therefore, 
it is believed that deep continuum surveys at >10 GHz 
may provide an excellent measure for the star formation 
activity in the highest redshift galaxies. 

This conclusion will only be valid if there is a negligi- 
ble amount of microwave emission associated with dust. 
While thermal emission from interstellar dust grains is 
negligible between 10 — 100 GHz, so called 'anomalous' 
dust emission, possibly arising from a population of ul- 
trasmall, rapidly rotating grains with a non-zero elec- 
tric dipole moment, may c ontribute significantly at m i- 
crowave wavelengths (e.g. iDraine fc L a zarian Ill998al fa: 
Ali-Haimoud, Hirata, & Dickinson 2009). If significant, 
high-frequency radio observations would likely overesti- 
mate true star-formation rates. This occurrence would 
also require that the observed high- z galaxies have a sig- 
nificant amount of ultrasmall dust grains. 

Comparing these observations with results from 
ALMA, which should be able to detect the peak FIR con- 
tinuum em ission from galaxies a t similar high redshifts 
(i.e. z >5: iMelchior et al. 1 1200 lh . should help disentan- 
gle this issue and result in powerful star formation and 
AGN diagnostics. Such tools will be critical for aiding in 
the interpretation of deep SKA imaging surveys which, 
unlike ALMA, will be able to survey large areas of the 
sky relatively quickly. And, as previously stated, even if 
the dust content in typical galaxies at these early epochs 
is negligible such that their FIR continuum emission is 
invisible to ALMA, radio continuum emission will likely 
be the best means to identify such objects and and quan- 
tify their star formation properties. 

4.4. Other Physical Considerations 

In the above presentation I have only considered vari- 
ations to the observed IR/radio ratios which arise from 
an increase in the competition of energy-loss processes 
to primary CR electrons relative to synchrotron radia- 
tion. Alternatively, other physical processes may have 
their hand in driving systematic variations in a galaxy's 
IR/radio ratio, perhaps even keeping the correlation lin- 
ear versus redshift fortuitously. 



4.4.1. Implications for the Possible Size Evolution of 
Galaxies 

It is possible that the suppression of a galaxy's non- 
thermal emission may not be as significant as previously 
suggested if the size of galaxies decrease with increasing 
redshift, resulting in lar ger values of Ur per u nit star 
formation. For instance. IBouwens et al. f (|2004h suggest 
that the half-light radius of galaxies between 2 < z < 6 
follow an rhi ~ (1 + z) -1 relationship. If true, the 
C^CMB oc (1 + z) 4 quenching of galaxy synchrotron emis- 
sion is effectively reduced by a factor of (1 + z) 2 for a con- 
stant magnetic field strength, and IC losses off the CMB 
will not be as significant for CR electrons. In this case 
such a strong departure from the canonical FIR-radio 
correlation may not be as easily seen. However, if galax- 
ies do become more compact, the escape times for CR 
electrons to leave the galaxy will likely decrease resulting 
in a decrease to the observed synchrotron emission. And, 
if the ISM of such galaxies are typically clumpy, similar to 
local dwarf irregular gal axies, then escape tim es should 
decrease even more (e.g. iMurphv et aT~l^00l) . 
For a random walk diffusion, the escape time is roughly 



proportional to l£ sc . Assuming that l c 



the escape time will decrease as (1 + z)~ 2 which should 
compensate for the increase to Ub by that same factor. 
Furthermore, if these compact galaxies also host star- 
bursts, advection out of each system due to a galactic 
wind may lead to even shorter escape times resulting in 
even larger IR/radio ratios at these early epochs. There- 
fore, even if galaxies are smaller at high-z, evolution in 
the FIR-radio correlation with redshift is still expected to 
be observed such that IR/radio should be larger with in- 
creasing z. Given that local starbursting galaxies follow 
the FIR-radio correlation, the suggestion for increased 
losses due to advection by galactic winds in such high-z 
systems would imply that they are different than star- 
bursts observed locally. 

4.4.2. Early Chemical Enrichment and the First Dust 

Grains 

The FIR-radio correlation may also fail at high red- 
shifts for a variety of reasons arising from differences in a 
galaxy's dust emission properties. For instance, IR/radio 
ratios should depart from the nominal value if dust is not 
present, the size distribution of grains take on a different 
form, or if the FIR emitting properties of dust grains are 
significantly different than in the local Universe. Cur- 
rently, observations of quasars and luminous galaxies at 
high-z sh ow that these syste ms contain large masses of 
dust (e.g. iWang et al. II2008D . While it is unlikely that 
these dusty systems are the representative population at 
early epochs, it is worth noting that the dust contained 
in these very young systems does not appear to be ex- 
otic enough to move systems off of the FIR-radio corre- 
lation. For example, the z = 6.42 quasar (J1148+5251), 
which is plot t ed in Figure [6] using data taken from 
iBeelen et al. I (|2006l) . h as an estimated dust mass rang - 
ing between 2 x 10 s M Q (|Dwek. Galliano, fc Jones 1120071) 
and 7x 10 8 M Q (jBertoldi et al. Il200l . Dust -to-gas mass 
ratios are also found to be > 0.004 which is similar to, or 
exceeding the, ~0.005 value of the Galaxy 

Using this QSO as an example, IDraine I (|2009f ) states 
that SNe are likely responsible for producing both the 
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metals which compose grains, as well as supernovae- 
condense d dust to provi de surface areas for grain growth. 
However. iDraine 1 (|2009f ) also shows that the dust in these 
objects are not predominantly composed of SN dust, 
known to have different emitting properties than that 
of typical dust grains in normal star- forming galaxies. 
Rather, the bulk of the dust formed in these high-z galax- 
ies is attributed to grain growth occurring more rapidly 
than grain destruction in the ISM, resulting in large dust 
masses that are already in place in galaxies 400 Myr after 
the first episode of star formation. This estimate assumes 
that star formation within galaxies began at z = 10. 
Given that the dust-to-gas ratios in these high-z galaxies 
do not appear significantly different than the Milky Way 
value, and that the grain formation process (i.e. in situ) 
is similar to th at of normal star-forming galaxies in the 
local Universe (|Draine 1120090 . it seems reasonable to be- 
lieve that the FIR dust emission from these galaxies does 
not drive large variations in their IR/radio ratios relative 
to the expected variations in their radio continuum emis- 
sion. 

4.4.3. Effect of CR Electron Acceleration Efficiencies 

The efficiency in which CRs are accelerated in SNRs 
plays a significant role in determining the synchrotron 
emissivity from galaxies. Furthermore, the distributed 
re-acceleration of CR electrons arising from external 
interactions such as ram pressure may add significant 
flux to the observed r adio continuum from galaxies (e.g. 
IMurphv et al. ll2009cD . The efficiency with which the ini- 
tial energy of the SNe explosion is converted into kinetic 
energy depends on the relative importance of radiation 
and adiabatic losses. In systems having extremely dense 
ISM, the adiabatic phase of the supernovae expansion 
may be halted as it expands into the ambient medium, 
thus reducing the amount of energy lost to adiabatic ex- 
pansion. It is possible that this extra energy could be, at 
least partially, used in the acceleration of CRs yielding a 
larger amount of synchrotron emission per unit star for- 
mation than in normal galaxies. For the SMGs, whose 
IR/radio ratios are systematically below the FIR-radio 
correlation by ~ 2 — 3 (sec ^4.2)1 . this scenario does not 
appear to help explain their low rat ios. 

As an example, the modeling of iDorfi I (|1991l . 120001) 
shows that the total acceleration efficiency for CRs in- 
creases from ~0.15i?sN to ~0.25.Esn, where Esn = 
10 51 erg is the total explosion energy of the SNe, when 
the external ISM density is increased from ~1 cm~ 3 to 
~10 cm~ 3 . The ISM densities of starbursts are likely 
much larger, as with the case of the example where 
niSM ~ 10 4 cm~ 3 was assumed, thus the evolution of 
SNRs will likely be different and may work to increase 
the synchrotron emissivity in such galaxies through an 
increased efficiency in particle acceleration. However, un- 
less SMGs host starbursts that are sufficiently different 
than those in local starbursting galaxies, which do not 
appear to exhibit such variations in their IR/radio ra- 
tios, this scenario alone might not be sufficient to explain 
their low IR/radio ratios. At present, more simulation 
work on this problem is needed. 

4.4.4. Secondary Electrons and Positrons 

So far synchrotron radiation arising for primary CR 
electrons accelerated in SNRs has been considered. CR 



nuclei, which are accelerated in the same SNRs, will in- 
elastically scatter off of the interstellar gas yielding tt° 
particles that decay into secondary e s and 7-rays, and 
7T^~ particles which will decay into secondary e^'s and 
neutrinos. These secondary e^'s will undergo the same 
energy-loss processes described in §2.1 and §2.2. At 
~GeV energies the ratio of secondary/primary electrons 
in the Galaxy is roughl y 2:1, and decrease s to nearly 
1:3 at ~3 GeV energies (jPorter et al~l 120081 ). The large 
secondary fraction arises from the fact that there are 
many more CR nuclei than electrons. Consequently, dif- 
fusion models suggest that ~20-50% of the diffuse GHz 
synchrotron emission in the Galaxy arises from secon- 
daries. However, secondary particles are produced more 
diffusively than primary electrons, thus their synchrotron 
morphology should not show any coincidence with their 
sites of origin. Given that the FIR an d radio images 
of galaxies are remarkably similar (e.g. IMurphv et al. I 
120061 120081 ) argues that the secondaries likely contribute 
mostly near or around star-forming regions. 

In the case of a starburst, having a much larger ISM 
density, the cross-section for collisions between CR nu- 
clei and the interstellar gas is increased, thus increasing 
the number of e^'s for a fixed primary nuclei/electron ra- 
tio which may actually dominate the diffuse synchrotron 
emission. Determining the contribution of additional 
synchrotron emission arising by an increased population 
of secondary e^'s clearly goes beyond the scope of this 
paper, however the fact that this effect may play a role 
is noteworthy. The recently launched Fermi 7-ray ob- 
servatory may help to shed some light on this topic. 

4.4.5. Variations in the Initial Mass Function 

While star-forming galaxies are currently thought to 
be responsible for completely reion izing the intergalac- 
tic medium (IGM) by z ~ 6 (e.g. iBecker et al~l I200H 
iFan. Carilli. fc Keating I I2006T) . the ionizing flux arising 
from star for mation in Lyman break g alaxies at simi- 
lar redshifts (|Bouwens et al. 1(20071 12008[ ) appears to fall 
a factor of > 6 below the minimum value required to 
maintain an ionized IGM for a given clumping factor 
and escape fraction under the assu mption o f a Sal peter 
stellar initial mass function (IMF). ICharv I (12001 has 
shown that this discrepancy can be reconciled by flat- 
tening the stellar IMF to have a slope of ~ —1.7 if reion- 
ization occurred at z — 9. The idea of a top heavy IMF 
at these epochs does not seem completely inappropriate 
given that low metallicity environments will favor the 
production of more high-mass stars. In this case, one ex- 
pects there to be both more SNe, increasing the amount 
of non-thermal synchrotron emission, as well as an in- 
crease to the ionizing photon rate per baryon, leading to 
an increase in the thermal radio continuum emission. It 
is worth noting that only a few SNe are required to en- 
rich a galaxy from primordial values to a few tenths solar 
metallicity, thus most galaxies observed with ages larger 
than a dy namical timescale are not likely to be that m etal 
poor (e.g. lWiersma et al. Il2009t lYaiima et al. Il2009h . 

Using the results of lCharyl (|2008l ). the ionizing photon 
rate for a stellar IMF having a slope of —1.7 is a factor 
of wl.7 times larger than a that for a typical Salpeter 
IMF slope of —2.3 over a mass range of 1 — 200 Mq. 
Similarly, for a given star formation rate, such a change 
to the slope of the stellar IMF should result in an in- 
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crease t o the core collapse supernova rate by a factor of 
«2 fe.g. lMadau. della Valle. & Panagia |[l9M) . Thus, if 
the stellar IMF in these high rcdshift systems transitions 
to being top heavy, the results presented here should not 
be significantly affected as both the thermal and non- 
thermal radio continuum output should change by fac- 
tors of <2. In the case that there is a significant amount 
of dust in these high-z systems, and assuming that their 
non-thermal radio continuum emission is completely sup- 
pressed by the additional IC losses to CR electrons by the 
CMB, their IR/radio ratios should still remain similar for 
a varying IMF given that both the FIR and thermal radio 
emission are directly proportional to the ionizing photon 
rate. However, since the increase to the ionizing photon 
rate will translate into an increase to the thermal radio 
emission, careful analyses of radio source counts at these 
rcdshifts, along with comparison between rest-frame UV 
studies, may help to illuminate possible variations to the 
stellar IMF. 

5. SUMMARY 

In this article I have put together a predictive anal- 
ysis for the expected evolution of the FIR-radio corre- 
lation versus redshift arising from variations in the CR 
electron cooling time-scales as IC scattering off of the 
CMB becomes increasingly important. As a result, I 
have looked at how discrepancies between observed and 
predicted IR/radio ratios can be used to characterize 
the magnetic field properties of galaxies at high-z, where 
the application of other observational indicators are ex- 
tremely difficult. In doing so, I have also looked at the 
value of deep radio continuum surveys in studies of star 
formation at high-z, particularly in the context of the 
SKA. My findings can be summarized as follows: 

1. Deep radio continuum observations at frequencies 
>10 GHz using next generation facilities like the 
SKA will likely provide the most accurate mea- 
surements for the star formation rates of normal 
galaxies at high z. The non-thermal emission from 
such galaxies should be completely suppressed due 
to the increased IC scattering off of the CMB leav- 
ing only the thermal (free-free) emission detectable; 
this situation may be complicated if 'anomalous' 
microwave emission from s pinning dust grains (e.g. 
iDraine fc Lazarian I Il9 98a) in such systems is not 
negligible. 

2. For normal star- forming galaxies to remain on 
the local FIR-radio correlation at high red- 
shifts requires extraordinarily large magnetic field 
strengths to counter IC losses from the CMB. For 
example, the magnetic field of a z ~ 6 galaxy must 
be >500 ptG to obtain a nominal IR/radio ratio. 
Thus, galaxies which continue to lie on the FIR- 
radio correlation at high-z, such as the sample of 
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radio-quiet QSO's, most likely have radio output 
which is dominated by an AGN, or the physics at 
work is simply unknown. 

3. The average IR/radio ratio among a sample of 
z ~ 2.2 SMGs is consistent with that for a sam- 
ple of z <~ 4.4 SMGs. This is in strike contrast 
with the factor of ~3 increase expected due to 
the increase in energy losses to CR electrons from 
IC scattering off of the CMB between these two 
redshifts. While this comparison suffers from low 
number statistics, it suggests that either the star 
formation driven IR and/or radio continuum emit- 
ting properties of these two populations of SMGs 
are different, or that AGN dominate the radio 
output of both of them. Simply assuming that 
their magnetic field strengths are significantly large 
(i.e. > 300 /iG) is problematic since, in the pres- 
ence of such strong fields, CR electrons emitting 
at GHz frequencies will cool primarily by ioniza- 
tion and bremsstrahlung processes rather than by 
synchrotron radiation which should lead to large 
IR/radio ratios unless other physics is at work. 

4. To detect typical LIRGs (L m > 10 11 L Q ) at all 
redshifts will require nJy sensitivities at GHz fre- 
quencies, specifically a 5cr point-source sensitiv- 
ity of «20 nJy (i.e. ctrms ~ 4 nJy). Thus, for 
the SKA to achieve this sensitivity for a reason- 
ably long (300 hr) integration necessitates that 
A eS /T sya ~ 15000 m 2 K _1 . At this sensitivity, 
the SKA will be sensitive to all galaxies forming 
stars at ~25 M yr _1 . This includes a significant 
amount of sourc es included in the UV l umin osity 
function work of iBouwens et al. I (|2007l I2008D be- 
tween 4 < z < 9. 

5. Deviations in the IR/radio ratios of normal star- 
forming galaxies at high redshifts will likely provide 
one of the only ways to put constraints on the ex- 
istence and strength of magnetic fields in the early 
Universe. For instance, if one knows the radio spec- 
tral index, galaxies with magnetic fields <50 fiG 
should be easily identifiable once beyond a redshift 
of >4 as their IR/radio ratios should be a factor of 
>4 larger than the nominal value, on average. 
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